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Most small, general purpose computers do not have hardware 
for the manipulation of floating point numbers. They have, 
therefore, implemented software packages to provide floating 
point functions. Two approaches can be taken to solve this soft- 
ware problem. 

The first approach is to provide a number of subroutines, 
each of which performs a specific function, e.g. floating addition. 
This approach requires the user to pass inputs in a standard 
manner to the subroutine. Since most small computers have only 
two accumulators, only one operand (which requires two words to 
represent) can be passed in the accumulators. The address of the 
second operand must be provided along with the subroutine call. 
A typical calling sequence requires two load accumulator 
instructions, a subroutine call, and an address word. Furthermore, 
one operand is always destroyed and replaced by the result (usually 
the accumulator operand) . This means that an intermediate result 
cannot be tested without destroying an operand — a function often 
useful in computation, h lengthy program requiring these sub- 
routines cannot make efficient use of storage. Further, the 
manipulation of operands becomes a tedious, cumbersome job. 
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The second approach is to provide an interpreter in which all 
the floating point functions are imbedded and which can simulate 
floating point registers in storage. The interpreter is given 
control by a subroutine call. Once in control, it accesses 
succeeding memory words and "interprets" a sixteen bit word in a 
manner similar to the hardware. These instructions are not 
executed in the hardware sense but are instead interpreted to 
provide extended machine features. 

This second approach was adopted for the NOVA floating point 
package and is the subject of this manual. The Interpreter 
provides four floating point accumulators which can be addressed 
and manipulated in a manner similar to the actual hardware 
accumulators. Floating point instructions are syntactically 
similar to machine instructions and are assembled in a similar way. 
For example, *the instruction 

ADD# $2f f 2 # SNR 
which adds AC0 to AC2, skips the next instruction if the result 
is non-zero, and does not change AC0 or AC2, has a similar floating 
point version. This instruction 

FADD# ^,2,FSNR 

* This manual requires the reader to have a good understanding 
of both, "How to Use The NOVA" and the Assembler Manual. 
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adds floating accumulator (designated as FhC0) to PAC2 # skips 
the next instruction if the result in non-zero, and does not 
destroy FhC0 or FAC2. 

The Interpreter has additional desirable properties. It can 
be implemented in Read Only Storage (ROS), and it is completely 
reentrant. These properties, in addition to the floating point 
instruction set, are described in the body of this manual. 
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2 • :--:-;uctur c of the Ir.te-.rgrot-:r 

The Interpreter is non -destructive , i.e. no instruction is 
modified in any manner during execution. Furthermore, it does not 
store temporary information within itself. Instead, it uses a 
writable area which roust be provided by the user. These properties 
enable the Interpreter to be wired in ROS as a customer option. 

Within the Interpreter, subroutine linkage is via a push 
down list. This list is maintained in the writable area provided 
by the user. This property, coupled with the properties already 
mentioned, make the Interpreter reentrant. This simply means that 
if a second user routine requires the Interpreter, it may interrupt 
the present routine, perform its function using the Interpreter 
and return to the first routine without affecting the state of the 
Interpreter. One implication of this is that the Interpreter may 
"call itself." This property has actually been used, enabling an 
extended version of the package to compute elementary trans- 
cendental functions using floating point instructions. 

The Interpreter requires eight page locations which the 

user must not destroy. These locations and their uses are; 

004 : will contain the starting address of the Interpreter 
005 z will contain the starting address of the initialization 
routine for the Interpreter 

006s one temporary word used by the Interpreter and saved by any 
routine reentering the package . 
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Wl i a word set up by the user to contain the address of a writable- 
area for the Interpreter 
040 z a word containing the address of a user written "get character" 

routine 
041s a word containing the address of a user written "put character" 

routine 
042-043; words containing linkage addresses to determine the 
location of extended code (if present) 
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3. Floating Point Number Representation^ 

Floating points numbers are internally stored in two 
consecutive 16-bit words. The form is; 




"S" is the sign of the mantissa, "M" , in bits 8-31. The mantissa 
is considered to be a normalized six digit hexadecimal fraction, 
and the range of the magnitude of the mantissa is therefore, 

16**-1 4 M-4 (l-16**-6). 
The characteristic, "C" , is the integer exponent of 16 in excess 
64 1$ code - The total range of magnitudes is thisj 

16**-1*16**-63^F^(1-16**-6)*16**63 
or approximately 
8.7*10**-7 8^F^7.2*10**7 5 
Any operand having a zero mantissa is represented in true zero 
from, i.e., bits 0-31 are 0. Negative numbers are identical to 
their positive counterparts except "S"=l instead of 0. 

The maximum error of a normalized mantissa is less than 
16**-6. 



* 
See also "How to Use The NOVA," Appendix C. 
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^ * Interprete r :.' :c ere! Cer re^r c : : tructur: 
4. 1, General 

The Interpreter provides four floating point accumulators. 
They are numbered 0,1,2, and 3 similar to the hardware accumulators. 
The designation FACn will be used for floating accumulator n. 
Arithmetic is performed accumulator to accumulator as with fixed 
point instructions. Operands can be accessed and stored using 
memory reference instructions. Instructions which reference 
floating point operands in memory should provide an address which 
points to the first word of the two word operand. If indexing 
is specified, the hardware index register is used. For example, 

FLDA 1,4,2 
loads FAC1 from two consecutive words in memory whose first word 
address is 4+ contents (AC2) . Certain instructions manipulate 
hardware AC2 and AC3 and will be described in Section 4.5. No 
facility is provided for manipulating AC0 and AC1 with the 
floating point instruction set. 
4.2. The Writable Area 

The basic Floating Point Interpreter requires 64 (decimal) 
words of contiguous writable memory. If the Extended Interpreter 
is used, 110 (decimal) words must be provided. The Extended 
Interpreter and its added features are described in detail in 
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Section 5. Other than the above differences in the length of 
the writable area, this description applies to both versions. 

The address of the first word of the writable area provided 
by the user must be stored in location 007 of page before any 
Interpreter commands are executed. If a second routine (or 
third, etc.) may reenter the Interpreter, this routine must 
provide an address pointing to a different memory area before 
reentrance is made. 

A number of flags are stored in this writable area which may 
be examined by the user. To access these flags, the user must 
exit from the interpretive mode (as described in Section 4.4). 
If an index register is loaded from location 001 , the user may 
access these flag words by an instruction of the forms 

LDA 0,n,i 
where i is the index register (2 or 3), and n is a constant 
displacement as described below. 

n = 0% The first word contains the overflow/underflow flags. 
If the result of an operation is less than 16**-64, 
bit 15 of this word will be set to indicate underflow. 
If the result of an operation is greater than 
(l-16**-6)*16**+63, bit 14 of the word will be set to 
indicate overflow. Other conditions may set overflow/ 
underflow. These will be discussed where appropriate. 
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This word is initially cleared. Once a bit is set, it 
will remain set. It is the user's responsibility to re- 
set these bits. 

n = 1: After an input conversion (see Section 4.5,3.) this 
word will be zero if no conversion was performed 
(because of an input error) . Otherwise the word will 
be non-zero. 

n = 2s After an input conversion, this word will contain the 
7-bit ASCII character, right justified, that served as 
the break character in the input stream. 

n = 3s After an input conversion, this word will be if 
no decimal point (**.") was encountered in the input 
stream. If a point was seen, the word will be 1. 

The diagram below summarizes these flags. 

\N R D 
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A convenient means of accessing these flags is to define four 
symbolic equivalences such ass 

FLGS = l overflow/underflow flags 

C0K = 1 ^conversion OK 

BRKC = 2 fbreak character 

Fp TF = 3 ? floating point flag 

Now if AC2 contains the address of the writable area, any flag 

may be accessed by statements like the following; 

LDA 0, FLGS, 2 

LDA l,COK,2 

LDA 1, BRKC ,2 

LDA 1, FPTF, 2 

4.3. Initialization 



The Interpreter, or more correctly the writable area, 
must be initialized before floating point instructions are 
executed. This initialization should be given once for every 

writable area to be used. The command is simply 

FIJSiI 
which generates the instruction 

JSR (§5 
(note that the initialization routine address is in location 
005), Location 007 must point to the writable area. This 
command destroys AC3, but preserves all other accumulators and 
the carry. The initialization routine clears the overflow/ 
underflow flags and sets up linkage for the push down list. 
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In tiring -.siting the Interpreter . .ode 

To use the Floating Point Interpreter, it is necessary 
to distinguish between the processing of normal instructions 
(described in "How to Use The NOVA") and the processing of 
floating point instructions. Since the latter instructions are 
"interpreted" (not executed per se) , the Interpreter roust be 
given control before floating point code will be executed 
properly. Whenever the Interpreter is in control, this will 
be called the "interpretive mode." Otherwise, the machine 
will be referred to as in the "normal mode." To enter the 
Interpretive mode, the command is 

FETR 
which generates the instruction 

JSR (£4 
As noted earlier, location 004 will contain the starting address 
of the Interpreter. Once in the interpretive mode, only the 
floating point instruction set can be used. (normal instructions 
will cause unpredictable results). To return to normal processing, 

the command 

FEXT 

must be given. 
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FETR destroys AC 3, but the contents of all other accumulators 
and Carry will be saved. AC2 and AC 3 can be used for indexing 
in the interpretive mode. Certain instructions also enable 
modification of the contents of AC2 and AC 3. FEXT will restore 
AC0 and Carry to their state before entering the Interpreter. 
AC2 and AC3 will reflect any changes caused by floating point 
instructions that modified their contents. 

The normal sequence of code using the Interpreter would be 
of the form; 

FIN I j INITIALIZE INTERPRETER 

LDA j^,CNST 

STA 0,TEMP ; NORMAL INSTRUCTIONS 
FETR ?ENTER INTERPRETER 

I FLOATING POINT INSTRUCTIONS 

FEXT ? EXIT INTERPRETER 

LDA 1,C5 

4.5. Floating Point Instruction Set 

This section contains a description of the instruction 

set for the basic Floating Point Interpreter. These instructions 

obey the syntactic rules described in the Assembler Manual. 

For example, all floating ALC instructions require a source 

accumulator designation and a destination accumulator designation. 

All floating point instructions begin with an "F" to distinguish 

them from the normal instruction set. Appendix A summarizes 

the instructions and Appendix B gives their octal encoding. 
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4.5.1. Memory Reference Instructions 

The instruction 

FLDA n,adr 

causes FACn to be loaded with the two woid operand at adr,adr+l. 

The instruction 

FSTA n ( adr 

causes FACn to be stored in memory at adr,adr+l. 

The instruction 

FJMP adr 

causes control to be transfered to the floating point instruction 

at adr. 

The instruction 

FJSR adr 

causes control to be transfered to the instruction at adr and 

AC 3 to be set to the value of the current location counter +1. 

For example , a floating point subroutine can be executed by the 

following: 

;MAIN PROGRAM 



FLDA 0,LOC 
FJSR SUBR 
FSTA 0,RSLT 



LOAD FAC^ 

JUMP TO SUBROUTINE 

STORE RESULT 



SUBRg 



j SUBROUTINE 
FJMP 0,3 | RETURN 
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The instruction 

FFIX adr 
causes the floating point number at adr, adr+1 to be converted 
to a fixed point, double precision integer (truncated) at 
adr, adr+i. If the conversion results in an integer whose 
absolute value is greater than 2**24-1, the overflow flag will 
be set and 2**24-1 will be returned as the magnitude. Sote 
that while floating point numbers are represented in signed- 
raagnitude format, fixed point values will always be represented 
in two's complement notation. 

The instruction 

FFLO adr 
causes the fixed point, double precision integer at adr, adr+1 
to be converted to a floating point number at adr, adr+1. 
Negative integers must be represented in two's complement format, 

The instruction 

FISZ adr 
causes the contents of adr to be incremented by one and the next 
floating point instruction in sequence to be skipped if the 
result is zero. 
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The instruction 

FDSZ adr 
causes the contents of adr to be decremented by one and the 
next floating point instruction in sequence to be skipped if 
the result in zero. FISZ and FDSZ should be used with fixed 
point, single precision integers — not with floating point 
numbers. 

The instruction 

FST3 adr 
causes AC3 to be stored at adr. 

The instruction 

FLD3 adr 
causes AC3 to be loaded from the contents of adr. FST3 and 
FLD3 operate on real accumulator 3. FST3 is useful for saving 
the return address inside a floating point subroutine. Note 
that the return address should always be saved if a floating 
point subroutine exits and then enters the interpretive mode, 
since FETR destroys AC3. These two instructions also provide a 
means for initializing loop counts without leaving the 
interpretive mode. 
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For example, 



FEXT 

LDA #,CNT 

STA 0,TEMP 
LP1; FETR 



FEXT 

DSZ TEMP 

JMP LP 1 

FETR 



can be replaced by 



LP1 



FLD3 CNT 
FST3 TEMP 



FDSZ TEMP 
PJMP LP1 



4.5.2. Arithmetic Instructions (ALC) 
4.5.2.1. Options 

The floating point ALC instructions are similar to 
normal ALC instructions. Two floating accumulators must be 
specified. The first is the source accumulator, the second 
the destination accumulator. 

Seven skip conditions are defined (in addition to the 
default "no skip"). These conditions are listed in Table 4-1. 
The conditions FSZR and FSNR should be used with caution. Since 
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floating point arithmetic is inherently approximate, the 
probability of obtaining true zero is very low. The normal 
procedure is to test a result (or difference) against a small 
quantity, £. For example, if we wish to test for the convergence 

of an iterative procedure we might use the following; 

EPSLN: 3544^ jEPSILON IS 2*16**-4 





FLDA 


0,ORSLT 


?GET OLE RESULT 


FLDA 


1,NRSLT 


?GET NEW RESULT 


FSUB 


1,0 


? OLD -NEW 


FLDA 


l f EPSLN 


?EPSILON 


FPOS 


0,0 


?ABS (OLD-KEW) 


FSUB 


1 # 0,FSLE 


rSKIP IF < EPSILON 



Tacl- £-1 



&AJ-J 111'; 5! T. o n I c 



Effect 



P8LT 
FSLE 

FSGT 
FSGE 
FSKR 
FSZR 
FSKP 



skip if result < 

skip if result <_ 0, 

skip if result > 0, 

skip if result >_ 0, 

skip if result f 0, 

skip If result - 0, 
unconditional skip 
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All ALC instructions permit the load/no load option. As 
with normal instructions, if a floating point instruction 
irnemonicis suffixed with "#", the results of the operation 
will not replace the contents of the destination register. 

A further option is available with one class of ALC 
instructions. This option will prevent post-normalization of 
the result. The instructions described in Section 4.5.2.2. 
permit this option. It is called for by suffixing "U" (for 
unnormalized) to the instruction mnemonic. For example, 
FMOV 0,1 moves FAC0 to FAC1 and normalizes the result, while 
FMOVU 0,1 moves FAC0 to FAC1 without normalization. 
4 -5.2.2. ALC Instructions with Post-Normalize Option 
The instruction 

FMOV n,m 
moves FACn to FACm. 

The instruction 

FPOS n # m 
moves the absolute value of FACn to FACm. 
The instruction 

FMNS n,m 
moves the negative of the absolute value of FACn to FACm. 



-19- 

The instruction 

FNEG n,m 
moves the negative value of FACn to FACnu 
The instruction 

FRND n,m 
rounds the value of FACn and moves it to FACm. By "round" we 
mean the following. 

F = i6**-6* L.16**6*F+1/2J 
The instruction 

FADD n,m 
adds FACn to FACm and moves the result to FACm. If underflow 
occurs, the underflow flag is set and true returned as the 
result. If overflow occurs, the overflow flag is set and a 
magnitude of ( l-16**-6) *16**+63 is returned as the result. The 
operands are assumed to be pre -normalized. 
The instruction 

FSUB n,m 
subtracts FACn from FACm and moves the result to FACm. The 
overflow conditions are handled as with FADD. Prenormalized 
operands are assumed. 

Table 4-2 summarizes the floating ALC instructions 
with post-normalize option* 



gives the maximum integer K such that K <^ X 



Table 4-2 



ALC Instructions with Post-Normalize Option 



Instruction 



Effect 



FMOV n,m 

FPOS n,m 

FMNS n , m 

FNEG n , ro 

FRND n,m 

FADD n , ro 

FSUB n,m 



FACn 



FACn 



-?- FACm 



-5*- FACm 



— (FACnl 



- FACn 



FACm 



FACm 



rounded FACn — ^ FACro 
FACn+FACm — ^ FACm 
FACm-FACn — 3>- FACm 
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4.5.2.3. ALC Instructions that Always Post-Normalize 

This class of ALC instructions always post-normalizes 
the result. They assume pre -normalized operands. Overflow is 
checked and indicated by setting the overflow flag and returning 
a magnitude of ( l-16**-6) *16**63 as the magnitude of the result. 
Underflow is checked and indicated by setting the underflow flag 
and returning true as the result. 

The instruction 

FMPY n,m 
multiplies FACn by FACm and moves the result Sf FACm. 

The instruction 

FDIV n,m 
divides FACm by FACn and moves the result to FACm. 

The instruction 

FHLV n f m 
halves FACn and moves the result to FACm. 

Table 4-3 summarizes the ALC instructions that 
always post-normalize the result. 



Table 4-3 



ALC Instructions that Always Post-rormalizc 



Instruction 



FMPY n # ro 
FDIV n,m 
FHAV n,ra 



Effect 



FACm*FACn -> FACm 
FACnn/FACn — > FACm 
FACn/2. > FACra 
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4.5.2.4. Floating ALC Instruction Examples 

Some examples of legal floating point ALC instructions 

are : 

PMPY 1,0 

FADD 0, 1,FSGE 

FSUB# 1,0,FSLT 

FMOVU 3 , 

FMOV 0,0,FLST 

FNEG 0,0,FSKP 

4.5.3. Input/Output Instructions 

The use of I/O instructions requires the user to provide 
two special routines. The first is an input routine which, when 
called, must return an ASCII input character, right justified in 
AC0 with bit 8 = 0. The address of this routine must be stored 
by the user in location 040 of page 0, 

The second routine is an output routine which, when 
called, must accept an ASCII output character, right justified 
in AC0 with bit 8 = 0. The address of this routine must be 
stored in location 041 of page 0. All output messages to this 
routine will be terminated by a null (all zero) character. 

These I/O routines must be reentrant for the Interpreter 
to be reentrant. (If the routines which interrupt and use the 
Interpreter do not use I/O instructions, the user I/O routines 
need not be reentrant) . 
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The instruction 

FDFC n ; FLOATING POINT DECIMAL TO FLOATING 
; CONVERT 

will cause an ASCII character string in engineering notation to 

be converted to internal floating point form and loaded in 

FACn. The input characters must be provided by the user routine 

whose address is stored in location 40 Q of page 0, Numbers in 

the following form will be converted. 



< l I- ; r ; L • • • ^ L-3 -n • — ?i[e W ^[^]j" 




where n is the decimal mantissa (the first seven non-zero digits 

will be converted and the remaining digits ignored) . The signs 

of the mantissa and characteristic are optional with the default 

assumed +. The break character is any character other than 

1. a decimal digit 
or 2. an "E" 

If the break character is a rubout (177), the entire 

string will be ignored and a new one must be given, i.e. the 

conversion starts over. 
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If the conversion results in a number less than 
16**-l*16**-63, the underflow flag will be set and true zero 
will replace FACn. If the conversion results in a number whose 
magnitude is greater than ( l-16**-6) *18**63 # this latter 
magnitude will replace FACn and overflow will be set. As 
described in Section 4.2, input conversion returns three 
additional words of informations conversion OK flag, the break 
character, and decimal point seen flag. These may be examined 
and used as necessary. 

Examples of legal character strings are; 

1* 
1.* 

-1* 

+1* 
IE 3* 

3.1415926* 
l.E+10* 
l.E-70* 

where * will be returned as the break character and conversion 
OK will be non-zero. 

Some illegal character strings ares 

A (break character will be A) 

+* (break character will be *) 
+.1 (break character will be I) 

Conversion OK will be zero in all these illegal cases. 
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The instruction 

FDFCI n ? FDFC WITH INDICATION 
will provide the user with an indication before the conversion 
begins. The ASCII character "F" followed by a null character 
will be passed to the output routine whose address is given 
in location 041. For example, if the user has provided for 
I/O from the teletype, the use of FDFCI will print "F" on the 
page copy every time an input is required. In all other 
respects it is identical to FDFC. 

The instruction 
FFDC n ; FLOATING POINT FLOATING TO DECIMAL CONVERT 
will convert the number in FACn to an ASCII character string 
in engineering notation. The output characters will be passed 
one at a time, right justified in AC0, to a user routine whose 
address is stored in location 41 Q of page 0. The output string 
will be of the form; 

( t I • T i ^ ^ w ^ ^ "H L |t 5 / m ^ 



where the n's represent the decimal digits of the mantissa, and 
the m»s represent the decimal characteristic. The string will 

be terminated by a null character. 
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4.5 S 4. special Instructions 

Two special instructions ari defined which modify the 
index registers. 

The instruction 

FIC2 
causes AC 2 to be incremented by two 

The instruction 

FIC3 
causes AC3 to be incremented by two. These instructions are 
useful for indexing through a table of floating point numbers 
(see Section 6, example #2) . 

A third special instruction provides a HALT feature within 
the interpretive mode, The instruction 

FHLT 
will cause the Interpreter to HALT. Hardware AC0 will contain 
the address of the FHLT instruction. The address lights will 
have no apparent relationship to the HALT, since the address 
is within the Interpreter. The user may press CONTINUE to 
resume after this HALT. 
4.5.5. Illegal Instructions 

The proper encoding for all floating point instructions 
is given in Appendix B. The Interpreter will HALT if an illegal 
instruction is encountered. Hardware AC0 will contain the 
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addr-~3 where the illegal instruction was found. This HALT will 
occur if extended instructions are used and only the Basic 
Interpreter is loaded , or on an} bit configuration that cannot 
be decoded into a floating point instruction. The user cannot 
press CONTINUE to resume after this HALT, 
4»6 Requirements for Reentrance 

A number of points regarding reentrance of the Interpreter 
have been mentioned. This Section explicitly defines the rules 
which must be followed by any routine which interrupts a base 
level routine and reenters the Interpreter. 

A. All hardware accumulators, Carry, and page locations 

006 and 001 must be saved, 

B. A new writable area address must be provided in 
location 001. 

C. An FINI must be issued after location 001 has been 
set up (only necessary the first time) . 

D» If I/O instructions are to be used, the user I/O 

routines must be reentrant . (Alternatively, locations 
040 and 041 must be saved and addresses provided to 
different I/O routines) . 

E« Upon exit to the base level routine, the hardware 
accumulators, Carry and locations 006 and 007 

must be restored. 
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5. .i:-:t c-ndef. ^lo^ting Foirt rcaturvS 

The instructions described up to this point are available 

with a IK Interpreter. An extended version of the Interpreter 

# 
is available which occupies 2K of storage. The extended version 

provides the instructions already described, in addition to a 

number of mathematical functions and ,! F" format output, If the 

extended Interpreter is used, 110 (decimal) words of writable 

storage must be provided by the user. 

5.1. Mathematical Functions 

The math functions are implemented using ALC instructions 
which always post-normalize (see Section 4.5.2.3.). They 
permit the no load option as well as the floating skip options. 
Appendix C provides a detailed description of the methods used 
to implement these functions as well as a discussion of their 
accuracy. The following is a general description of each 
instruction. 

The instruction 

FALG n,m 

computes the natural logarithm of the contents of FACn and 

moves the result to FACm. If the argument is less than f 

the overflow flag is set and -( l-16**-6) *16**63 is returned as 

the result. 

With a 4K configuration, the Interpreter requires locations 
5600-6577 (octal) and the Extended Interpreter requires 
locations 4100-5577. 
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The instruction 

FATN 11 , m 
computes the arctangent of the contents of FACn and moves the 
result to FACra. The result is an angle expressed in radians 
in the range -pi/2 <C arc tan (x) ^ pi/2. 

The instruction 

FCOS n r m 
computes the cosine of the contents of FACn and moves the result 
to FACm. The argument is assumed to be an angle expressed in 
radians. If the argument is greater than 2**24, the overflow 
flag is set and the result will be incorrect. 

The instruction 

FSIN n,m 
computes the sine of the contents of FACn and moves the result 
to FACm. The argument is assumed to be an angle expressed in 
radians. If the argument is greater than 2**24, the overflow 
flag is set and the result will be incorrect. 

The instruction 

FTAN n,m 
computes the tangent of the contents of FACn and moves the 
result to FACra. The argument is assumed to be an angle 
expressed in radians. If the argument is greater than 2**24, the 
overflow flag will be set and the result will be incorrect. 
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The instruction 

FEXP 11 »ro 
computes e raised to the power contained in FACn and moves 
the result to FACm. If the argument is less than -177.5, true 
is returned and the underflow flag is set. If the argument 
is greater than 174.673, +( l-16**-6) *16**63 is returned and the 
overflow flag is set. 
The instruction 

FSQR n,m 
computes the square root of the argument in FACn and moves 
the result to FACm. If the argument is less than 0, the under- 
flow flag is set and true is returned as the result. 

Table 5-1 summarizes the math functions. 
5.2 . "F" Format Conversion 

The standard Interpreter provides floating point to decimal 
conversion with "E" format output. The extended version provides 
"F" format output as well. 
The instruction 

FFDCF n r FFDC WITH "F" FORMAT OUTPUT 
will convert the floating point number in FACn to decimal and 
output a character string in "F" format via the user routine 
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whose address is stored in location 041. The output will be 
of the form: 



a a a \ f\% 



^H* ^71'Ti 1\ 



w 



The width of the field (including sign and decimal point), W, 
and the number of places to be given after the decimal point, 
D, roust be set up in the writable area before FFDCF is given. 

The displacements of these words are 

W = 136 
D = 137. 

They can be accessed in a manner similar to the flags described 

in Section 4,2. Two conditions will cause the overflow flag 

to be set and no conversion to be performed. 



Table 5-1 



MATH FUNCTIONS 



Instruction 



FALG n,m 

FATN n , m 

FCOS n,m 

FSIN n # m 

FTAN n , m 

FEXP n , m 

FSQR n,m 



Effect 






In(FACn) 


^ 


FACm 


arctan (FACn) - 




-FACro 


cosine (FACn) - 




-FACm 


c i r%a ^ PSTn^ 




- FACm 


bine \Er\\~llf 

tangent (FACn) 




-FACm 


e** (FACn) 


— >. 


.FACm 


(FACn)**^ — - 




FACm 
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1.) K >32 (entire width of field limited to 32 characters) 

2.) w < D+2 (w must be 2 greater than D to provide for 
sign and decimal point) 

If W is not large enough to accommodate the number, significant 

digits will be lost. 

Assume W has been set to 12 (decimal) and D to 6. Examples 

of "E" format versus "F" format outputs are; 



L^_Fprrr,at,. 



+.1374600E+02 
-.7968433E-03 

+.1000000E+04 
-.1000000E+05 
+.1000000E+06 
-.3500000E-06 
+.4713279E-01 



F" Format 



+13.746000 

-.0007 96 
+1000.000000 

10000.000000 (note sign lost) 
00000.000000 (all lost) 

-.000000 (significance 
+.047132 lost) 
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fl ~ Examples. 

The following ,vont-JL»e is an example of a floating point 
subrom-ine that performs a square root Newton Iteration on the 
trial guess in FUCGf given the argument in FAC1. 



THIS ROUTINE PERFORMS A SQUARE ROOT 

NEWTON ITERATION 
FAC1 CONTAINS THE ARGUMENT AND IS 

NOT DESTROYED 
FAC0" CONTAINS THE PREVIOUS GUESSTIMATE 
COMPUTES (FAC0+(FACl/FACJ2T))/2. 
CALLING SEQUENCE 

FJSR NSR 

RETURN : RESULT IN FAC0 



NSR: FMOV 1,2 ;SAVE ARGUMENT 

FDIV 0,2 ?FAC1/FAC0 

FADD 2,0 ,-FAC0+FACl/FAC^ 

FHLV 0,0 ?(FAC£f+FACl/FAC£f)/2. 

FJMP 0,3 ; RETURN 

The next routine is an example of polynomial evaluation. 
This routine requires AC2 to point to the first word of a table 
of floating point coefficients, ordered high order coefficient 
down and terminated by true 0. For example: 

AC 2 r~~ ^^ 



CJX 



Cn-1 



Cn-2 



Word 


1 
2 
3 

4 



-, 
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Cl 



Co 







JL 







SL 



Word 



2n+^ 
2n+l 
2n+2 
2n+3 



The routine uses Horner's method for evaluation, i.e. 
f(x) = (•••(((x + C n ) x + c n-1 ) x + C n _ 2 ) x---+ Cj) x + C Q 
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t OT.YtlOMTUT, EVA I,UAT ION 

FAC2 CONTAINS ARGUMENT (X) 

AC 2 POINTS TO COEFFICIENT I,T/7T TERMINATED BY 0, 

AND ORDERED HIGH TO T.OW 
RESULT RETURNED IN FAC0 
FAC0,FAC1 DESTROYED 
CALLING SEQUENCE 

FJSR FPLY 

RETURN 



fCLEAR RESULT 
?GET COEFFICIENT 



FZRO: 









FPLY: FLDA 


0,F2RO 


FPLY1: FLDA 


1,0,2, 


FMOV 


1,1,FSNR 


FJMP 


0,3 


FMPY 


2,0 


FADD 


1,0 


FIC2 




FJMP 


FPLY1 



; RETURN IF ZERO 

I SUM * ARGUMENT 

I SUM * ARG. + COEF. 

; BUMP POINTER TO NEXT COEF 



APPENDICES 



Appendix £ 



Flo: ting P oint Instru-t icn r ur.T.?r; : 

Standard 



FftDD 

FDFC 

FDFCI 

FD1V 

FDSZ 

FETR 

FEXT 

FFDC 

FFIX 

FFLO 

FHLT 

FHLV 

FIC2 

FIC3 

FINI 

FISZ 

FJMP 

FJSR 

FLD3 

FLDA 

FMOV 

FKNS 

FMPY 

FNEG 

FPOS 

FRND 

FST3 

FSTA 

FSUB 



Floating £dd 

Floating Decimal to Floating Ccnvart 

FDFC with"~_Indication 

Floating Divide 

Floating Decrement and £kip if Zero 

Floating Mode Enter 

Floating Mode Exit 

Floating Floating to Decimal Convert 

Floating to Fix ed 

Fixed to Floating 

Floating Halt 

Floating Halva 

Floating Increment AC2 

Floating Increment AC3 

Floating Initialize 

Floating Increment and Skip if Zero 

Floating Jump 

Floating Jump to Subroutine 

Floating Load £C3> 

Floating Load Floating Accumulator 

Floating Mov e 

Floating Mov3 Minus 

Floating Multiply 

Floating Negate 

Floating Move Pos itive 

Floating Round 

Floating Store AC3 

Floating Store Floating Accumulator 

Floating Subtract 



5xt:::nd.jd 

FBLG 
FATN 

FCOS 
FEXP 
FFDCF 

FSIN 
FSQR 
FTAN 



Floating Natural Logarithm 

Floating Arctangent 

Floating Cos ine 

Floating Exp onential 

Floating Floating to Decimal Convert 

~~ with "F" Format 

Floating Sin e 

Floating Square Root 

Floating Tangent 



Floating Point Options 



# 
FSGE 
FSGT 

FSKP 
FSLE 
FSLT 
FSNR 
FSZR 
U 



No Load 

Floating Skip on Greater Than or Equal 

Floating .Skip on Greater Than 

Floating Skip 

Floating Skip on Less Than or Equal 

Floating £kip on Less Than 

Floating jSkip on Non Zero Result 

Floating Skip on Zero Result 

Unnorraalize (no post-normalization) 



Appendix B 
"_"."■■ -. \r-. -r "■■.:.!r- >;truai:r Prcc-i-a 
i'.enory Reference 

J 2.3 4.:H7 9 



IS 



— f— — 
'A r- 4 


i 
1 1 

1 1 


h — P_ 

f 
1 


D 



Fl 

Fl 



01 

10 



LDA 
STA 



F2 is FAC 



Fl = 00 



F2 


= 


00 


FJMP 


F2 


= 


01 


FJSR 


F2 


= 


10 


FISZ 


F2 


= 


11 


FDSZ 


F2 


— 


00 


FFLO 


F2 


= 


01 


FLD3 


F2 


=s 


10 


FST3 


F2 


= 


11 


FFIX 



Fl 



11 



Instructions Requiring an Accumulator 



12 3 4 5 



N 15" 




Fl = 01 

Fl = 10 
Fl = 11 



FDFC F2 = 1 "I" 
FFDC F2 = 1 "P" 
Illegal 



ilal Instructions 









3 4 


s 










15 






100 |f t|0- 




_ .. - __ rf\ 






■ i 


Fl = 00 
Fl = 01 
Fl = 10 
Fl = 11 








FEXT 
FIC2 
FIC3 
FHLT 


7. LI Instruction:; 










y; i ; .? 4 s 7 


? 


9 m w \i 


13 15 






1 

* 

i ACS 
_J , 


ACD 


fi 


— 

F 


1 


F3 * 

[ ._ 


I 


S = 

S = 1 
S = 2 
S = 3 


No skip 
FSGT 
FSLT 
FSNR 


► 










S 

s 
s 

F 


ss 


4 FSZR 

5 FSGS 

6 FSLE 

7 FSKP 



Fl ^ (8-10 must be 0? 11 is "u" option) 



Fl 


rs 


1 


FNSG 


Fl 


= 


2 


FMOV 


Fl 


= 


3 


FPOS 


Fl 


= 


4 


FMNS 


Fl 


= 


5 


FSUB 


Fl 


= 


6 


FADD 


Fl 


= 


7 


FRND 



= 










F3 = 










F2 = J 


. 


FMPY 




F2 = 2 


> 


FDIV 




F2 = 2 


) 


FHL¥ 




F3 ]4 (Extended) 




F2,F3 


= 


1 


FALG 




F2 # F3 


= 


2 


FATN 




F2,F3 


= 


3 


FCOS 




F2, F3 


ss 


4 


Illegal 


(FMPY) 


F2, F3 


= 


5 


FSIN 




F2 # F3 


= 


6 


FTAN 




F2,F3 


= 


7 


Illegal 




F2,F3 


= 


10 


Illegal 


(FDIV) 


F2 # F3 


= 


11 


FEXP 




F2,F3 


= 


12 


FSCR 




F2,F3 


ss 


13 


Illegal 




F2,F3 


= 


14 


Illegal 


(FHLV) 


F2 # F3 


ss 


15 


Illegal 




F2,F3 


ss 


16 


Illegal 




F2,F3 


= 


17 


Illegal 





Appendix C 

The following pages describe in some detail the function, 
method, and accuracy of the mathematical routines supplied 
with the extended Interpreter. The accuracy of the routine 
is influenced by two factors: 

1. The accuracy of the argument and 2. The accuracy of 
the algorithm. These two factors will be mentioned for each 
routine. The accuracy of the algorithm itself assumes an 
argument that is exact, i.e. no argument error. 

The relative and absolute errors of a function routine 
are defined as follows: 

Let f(x) = the true value of the function at x 

Let g(x) = the result returned by the function routine 
given x 

Now the absolute error of the result is, 

ABS(f(x) - g(x)) 
and the relative error of the result is, 

ABS((f(x) -g(x))/f(x)). 



Floating Point Arctangent 
Functions To calculate the arctangent of x, where x is a 
floating point number, and return an angle in 
radians in the range -pi/2 ^a re tan (x)^pi/2 
Method; The range of x is immediately reduced to 

by means of the identities 

arctan (-abs(x)) = - arctan( abs(x) ) 
and if abs(x) > l r 

arctan abs(x)) = pi/2 -arctan ( l/abs(x) ) 
For x > tan (pi/12) the range is reduced to 

tan (-pi/12) <. y ^ tan (pi/12) 
by means of the identity 

arctan (x) = arctan (x0) + arctan (x-x0/( l-x*x#) ) . 
For x0 = 1/3**0.5, we obtain 

arctan (x) = pi/6 + arctan (x*3**0.5-l/(x+3**0. 5) ) 
= pi/6 + arctan (y) 

w here y satisfies the range given above. The arctan 
is computed using the first four terms of a poly- 
nomial approximation of the form: 

arctan (x) = X * -§E Ci*x**(2*i) 



Accuracy : 
Argument Error 

If x is the argument, the absolute error of the result 
is approximately 

£/(X+x**2) 
where £, is the absolute error in x. Thus for small values 
of x, the errors are almost equal, while as x becomes larger, 
the effect of the argument error decreases. 
Maximum Relative Error 

For the range 

-tan(pi/12) ^ x 4. tan(pi/12), 
the maximum relative error is approximately 10**-7.6. 



Floating Point Exponential 

Function: To calculate e to the power x, where x is a floating 

point number, 
Methods If x <-177.5 # return as the result and set underflow 

flag. 

If x> 17 4. 67 3, return the largest positive number as 
the result and set overflow flag. 
Otherwise, we use the equality 

e** x = 2**{x*log2(e) ) 
Let x*log2(e) = m*16**c 

Further, let q = Lm*16**cJ 

and f = m*16**c-q 
Therefore J2f^f <1.0 

Now 2**(f+q^ = 2**f*2**q 
Let us compute 2**f. 

The range of f can be further reduced if we let 
g = f - 0.5 where -0.5 ^g ^0.5 
We compute 2**g directly if g >0, otherwise we 

compute 2**g as 

l/(2**abs(g}} 

2**g (where < g <#. 5) is computed using the first 
five terms of a polynomial approximation of the forms 

ri 
2**g ~ 2E cx*g**i 



Now f = g+0.5 

Therefore, 2**f = 2**g*2**0.5 

But the answer is 2**f*2**q 

q is an integer and we let q = 4i + j 
Now 2**f*2**q = 2**f*16**i*2**j 

The characteristic of f is added to i to obtain 
2**f*16**i. This result is shifted j positions 
left if i>0 or j position right if i<0. 
The result is of course e**x. 
Accuracy; 
Argument Error s 

The relative error of the result is approximately equal 
to the absolute error of the argument. Thus for large values 
of x, substantial relative errors in the results can occur. 
Maximum Relative Error 
For the range 

0^x <0.5 
the maximum relative error of 2**x is approximately 10** -7.0. 



Floating feint l'.~tur?i Log^ri* hro 

Function; To calculate the natural logarithm of x # where x i?: 

a floating point number. 
Methods If x<^, overflow flag is set and minus the largest 

floating point number is returned as the value. 

Otherwise, let x = m*16**c. By means of a binary 

normalization, the range of m' is reduced to 
1/2 4m 1 <1, 

and x = m* *16**p*2** ( -q) where g = number of left 

shifts required to normalize (0 ^.q <*3) . 

Now for l/2( v m' <l/2**{2f # 5 

let a = 1/2, b = 1 

for 1/2^*0.5 ^ m 1 < 1 

let a = 1, b = 

Define y = (m* -a)/(m* +a) 

then m 1 = a*(y+l)/( -y+1) 

Now x = 16**p*2**(-q)*a*(l+y)/(l-y) 
= 2**(4p-q-b)*(l+y)/(l-y) 

Using ln(x) = ln(2) *log2(x) we obtain 

ln(x) = ln(2)* [(4p-q~b)+log2((l+y)/(l-y))3 

ln(x) = ln(2)* (4p-q-b)+ln((l+y)/(l-y)) 

From the above, we can determine that 
2**H2f*5 ^ (l+y)/(l-y)^2**0.5 



The ln( (l+y)/(l-y) ) is determined for the above 

range using the first three terms of a polynomial 

approximation of the for"; 

ln(z) = z* <>_ ci*z**(2*i) 

where z = (l+y)/(l-y). 
Accuracy; 
Argument Error 

The absolute error in the result is approximately 
equal to the relative error in the argument. Therefore, an 
argument close to 1 can give a large error since the function 
at this value is quite small. 

Maximum Ab s olute Error 

For the range 

l/2**-0.5 ^ (l+x)/(l-x) 4.2**0.5 

the maximum absolute error of ln(x) is approximately 
10** -7. 6. 



Floating Point c irr_ ~nd Cosine 



Function: To calculate sin(x) or cos(x), where x is the 

floating point angle in radians. 
Method; Compute p = abs(x)*4/pi 

Let q = )p|, f = p-q where 0<f<l 
Now q represents the half-quadrant in which the abs(x) falls. 
Using the following equalities* 

sin(x) = -sin(-x) 

cos(x) = cos(-x) 
cos(x) = sin(x+pi/2) 

we define ql = q if sin is required and x%0 
ql = q +2 if cos is required 
ql = q + 4 if sin is required and x <0 

Then for all values of x, the computation has been reduced to 
sin(pi/4*(q +f)) = sin(t) 

Since sine ( and cosine ) are periodic in 2*pi, we take qj_ = q^ rood 8 
Using the further equality 

sin(pi/4+x) = cosine (pi/4-x) 
we finally can produce the table below 
q x sin(t) 




1 
2 

3 

4 

5 
6 
7 



sin(pi/4*f) 
cos(pi/4*(l-f)) 
cos(pi/4*f) 
sin(pi/4*(l-f)) 

-sin(pi/4*f) 
~cos(pi/4*(l-f) } 
-cos(pi/4*f) 
-sin(pi/4*(l-f)} 



I~ ""ill •/?."- ~, the argument lias been reduced to the range 

0^,t <pi/4. The sin is computed using the first four terms of 

a polynomial approximation of the form; 

sin(x) ~ x* 2* Ci"*x**(2*i) 
The cosine is computed using the first four terras of a poly- 
nomial approximation of the forms 

n 

cos(x) ^ 2E: Ci*x**(2*i) 
i ~<Z) 
Accuracy; 

Argument Error 

The absolute error of the result is approximately 
equal to the absolute error in the argument. Thus, the larger 
the argument, the larger the absolute error of the result. 
M aximum Relative Error 

For the range 

^ x < pi/4 
the maximum relative error for sin(x) and cos(x) is approximately 
10** -7. 4 



ric^tir.g Point Square Root 
Function: To calculate the square root of a floating point 

number, x. 
Method: Let x = m*16**c 

If m<0, set underflow flag and return 0. as the 

result. 

If m=0, return h as the result, 

Otherwise, let c = 2*p + q where p is an integer 

and q = or 1. Now if q = 0, we have x = m*16**2p 

and x**l/2 = ro**l/2 * 16**p. 

If q = 1, we have x = ra*16** ( 2p+l) 

= ro*16**(2p+2)/16 

and x** 1/2 = (m**l/2)/4 *16**(p+l) 

Therefore, the characteristic of the result is 

p + q, and the problem has been reduced to finding 

a suitable first guess for nn**l/2 if q = or 

(m**l/2)/4 if q = 1, 

An initial guess is taken in the hyperbolic form 

Y0 = a+b/(c+ro) 



where for q = 0, 



and for q = 1, 



a = 1.80713 
b = -1.57727 
c = 0.954182 

a = 0.428795 
b = -0.3430368 
c = 0.877552 



The ir It ill gjnss is now 

yl = y^*16**(p+c;) 
Two foewton iterations give us the result. 

y2 = (yl+x/yl)/2 

x**l/2 = y3 = (y2 + x/y2)/2 

Accuracy; 

The relative error of the result is approximately 
half the relative error in the argument. 
:-^2£.i^lIl__^Ll_^i v? Error 

The maximum relative error for x**0„5 is approximately 



Floating I-'cint Tangent 

Funrt-iori! To calculate the tangent of x, where x is the floating 

point angle in radians. 
Methods Compute p = abs(x)*4/pi 

Let q = [pj , f = p-q where < £ <1 
Take q = q rood 4 and in a manner similar to sine- 
cosine we can obtain the table below 
q^_ tan(x) 




1 
2 
3 



tan (pi/4*f) 
cot (pi/4*(l-f)) 
-cot (pi/4*f) 
-tan (pi/4*(l-f)) 



In all cases, the argument has been reduced to the 
range 0"^arg<pi/4. 

The tangent is computed using the first six terras of 

a polynomial approximation of the form: 

n 
tan(x)= x* iS Ci*x**(2*i) 

< = 

Accuracy: 
Argumen t Error 

The absolute error of the result is approximately 
equal to 

& * (1 + tan(x)**2) 

where £ is the absolute error of the argument. Thus if x is 
near an odd multiple of pi/2, an argument error will produce 
a larger absolute error in the result. 



-' - -: :< 1 r 1 1: *n -» - j a t i ve Error 

For the range 

0^.x<pi/4 

the maximum relative error of tan(x) is approximately 10** -6. 6, 



